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ABSTRACT 



Two models for the correlation of ship position reports obtained 
and relayed by ocean surveillance satellites are described. The models 
are intended to process position reports of ships detected by satellites 
and correlate this data to provide intelligence concerning the position, 
course, and speed of ships traversing the ocean surface. A probabilis- 
tic model assigns position report correlations by reliance on probabil- 
ities assigned to computed ship speed and to changes of course and 
speed. A reliability factor is determined for each track. The deter- 
ministic model avoids explicit use of probability theory and correlates 
position reports on the basis of proximity of the position report to a 
position predicted by track extension. Evaluation of the models by means 
of computer simulations is described. Results of the evaluation are 
discussed and recommendations are made concerning areas for further 
investigation. The computer programs are included in the Appendicies. 
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1. Introduction 



Since the first sputnik was successfully projected into orbit in 
1957, the possibility of using a satellite for reconnaissance and 
surveillance of the oceans has been of increasing interest to military 
men. Studies have been conducted to analyze the usefulness and feasibil- 
ity of a surveillance satellite system and several reports have proposed 
system requirements. A report by the Pacific Missile Range described 
some of the problems to be encountered in processing data received from 
a surveillance satellite, (VJ 

During the spring and summer of 1964, analysts at the Advanced 
Programs Division of the Pacific Missile Range conducted a computer 
simulation study employing a surveillance satellite to detect and report 
the positions of ships at sea. The purpose of the study was to determine, 
for a given satellite sensor capability and orbit, the times and positions 
of ships at sea eligible for detection by the satellite. The ship data 
used in the study was extracted from actual weather reports. (The studies 
are reported in detail in Pacific Missile Range Technical Note No. 3285- 
570, 1964.) The Pacific Missile Range program permitted the user to 
describe the orbit of a satellite in altitude and inclination to the 
earth's equator. It also permits the analyst to specify several different 
sensor capabilities. Ship's call signs, reported position, time, and 
cloud cover were extracted from the weather reports, and from this, tracks 
were projected. The simulation program then determined the times and 
positions of the ship at which the satellite would be in a position to 
detect the ship with each of the sensors described by the user. The 
simulation program also determined whether this conjunction of ship and 
satellite occurred during daylight or darkness. A Monte Carlo process 
based on reported cloud cover was used to determine whether the ship 



1 



was obscured by clouds at the time of conjunction. 

The object of this thesis is to provide the next logical step in 
the analysis of a surveillance satellite system, i.e., to devise a 
correlation logic to determine to what extent the sets of conjunctions 
provided by the Pacific Missile Range simulation can be processed to 
form tracks which correspond with the actual track as determined from 
the ship's weather reports. It will then be possible to determine the 
coverage frequency required to provide an acceptable tracking reliabil- 
ity, hence the number of satellites required. The authors have there- 
fore addressed themselves to the task of devising correlation logic to 
provide programs which will facilitate further study of the feasibility 
and cost of tracking ships at sea by use of data collected by surveillance 
satellites. 
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2. General Analysis 

The formulation of track correlation logic which will facilitate 
testing the feasibility of satellite surveillance of ocean shipping, 
has been undertaken via two different approaches. The two procedures 
considered are described in detail in the following two chapters. In 
the development of both procedures it is assumed that the input data 
to the tracking program will consist only of the time of conjunction 
and the geographical position of the contact at the time of conjunction. 

Thus the analysis is not dependent upon the development of a high reso- 

.4*: ' 

lution sensor which might provide information about the ship's course 
by virtue of its orientation to the orbital path, or other identification 
clues, such as length or the number of stacks observed. 

This study has been limited to a consideration of methods by which 
the data described above might be correlated in a computer program to 
derive tracks of ships at sea. No attempt has been made to determine 
optimal orbits or combinations thereof. 

It is evident that ocean traffic intelligence gathered from other 
sources would be of significant value in decreasing the uncertainty 
involved in surveillance by satellites. Such other sources include 
sighting reports by aircraft and friendly shipping, reports of ship 
arrivals and departures gleaned from commercial sources or news media, 
and voluntary reporting of own position by ships at sea. If one half 
of the ships in a certain area can be identified by means of external 
intelligence, the number of possible correlations to be determined may 
be reduced to one quarter of what would otherwise be required. This 
report will not concern itself with the methodology or effect of such 
an input to the tracking program, for if tracking can be shown to be 
feasible using information gathered solely by the surveillance satellite, 
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it will clearly be enhanced by inclusion of externally gathered 
information. 

2.1 The Correlation Problem 

The correlation problem can be thought of as consisting of two 
parts. The first and most difficult part involves track initiation, 
i.e., making a decision as to which pair of positions reported on 
consecutive orbits form the most probable true correlation. Once this 
decision has been made, the track can be considered as true and the 
course and speed of this designated track then becomes available for 
the second part of the problem, the track extension. 

Track Initiation . The speed with which ships sail the surface of 
the oceans is less than some generally known maximum speed. Few ships 
today can sustain a speed in excess of 40 knots for more than a few 
hours. An analysis of merchant ship cruising speeds reveals that their 
distribution is approximately normal with a mean of 14.0 knots and a 
standard deviation of 2.5 knots. Only 2.0% of the cruising speeds ex- 
ceeded 19 knots while 78.5% were in the range 11.5-19.0 knots. The 
inclusion of military ships in the analysis presumably would bias the 
results toward somewhat higher speeds but even fast naval ships seldom 
cruise at speeds in excess of 20 knots. Thus it may be reasoned that 
a correlation between two position reports which indicates a speed of 
16 knots is much more likely to be a true correlation than one which 
indicates a speed of 10 knots or 25 knots. The existence of this speed 
limitation determines a selection rule for track initiation which per- 
mits the elimination of all possible correlations between position 
reports which would require a speed in excess of some specified maximum 
speed. In many cases this selection rule alone will determine the cor- 
relation. i.e., if for a point from one sweep there is only one 
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corresponding point from the next sweep which is within a circle whose 
radius is equal to the maximum speed attainable multiplied by the time 
between the position reports, then there is a high probability that 
these two points can be correlated as a true track. It is evident that 
any tracking program must not deny correlations showing very high or 
very low speeds, but much more confidence can be placed on a correlation 
which indicates a speed in the range of 12-18 knots. 

The difficulty of correlation increases when a number of possible 
correlations occur within this maximum distance circle. The number of 
tracks to be considered as possible for position reports from two orbits 
is equal to the number of points from the second orbit which lie within 
the maximum distance circles drawn about the points from the first orbit. 
This situation is illustrated in Figure 1. In areas of high shipping 
density it is to be expected that the picture will be much more compli- 
cated than that illustrated. 




Note: The notation p rn refers to report number n on orbit number r. 

Illustration of Possible Tracks for a Given Situation 

Figure 1 
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When the circles of maximum attainable distance are drawn in 



accordance with the selection rule for track initiation it can be 



seen that while p may link with all four points obtained on orbit 

two, the same is not true for p, , p, , and p . Thus it can be 
’ 12 ’ 13 ' 14 

assumed that since p can only correlate with p . then p cannot 

14 21 ’ *21 

be logically considered as a possible correlation for p . A similar 

reasoning applies to p , and p , If p is removed from considera- 

13 24 24 

tion with p. , then p, can be linked with p , The result of apply- 
12 ’ ^12 *23 

ing the selection rule on the position reports illustrated in Figure 1 
therefore reduces the eight possible tracks illustrated to the following 
firm correlations; 
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This selection rule cannot stand alone however, since it is 
heavily dependent upon the assumption of a sensor which is 'perfect' 
in the sense of a radar with a blip-scan ratio of one, i.e., if a 
ship is within the sweep path of the sensor it will be detected by the 
sensor with probability one. It cannot be strictly applied when the 
circle of maximum attainable distance drawn about a point encloses a 
seaport since there is a possibility that the ship did enter port and 
was not detectable by the satellite. Similarly, any ship detected within 
a circle of maximum attainable distance drawn about a seaport may be a 
ship just starting her voyage from that port. The rule also ignores the 
possibility of a temporarily surfaced submarine or occasional data points 
resulting from contacts with aircraft or heavy clouds. 
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Track Extension , After tracks have been initiated, a rule of 
consistency utilizing both course and speed can be applied to determine 
the extension of tracks for position reports from consecutive orbits. 

This rule is one that states the extension of initiated tracks to new 
position reports can be done on the basis of course and speed consis- 
tency. Ships at sea are generally steaming from one port to another 
port and have selected a direct course. This implies that most course 
changes in general are executed to circumnavigate points of land or are 
occasioned by circumstances such as avoidance of severe weather, mid- 
course changes in destination, or navigational course adjustments 
necessitated by the effects of wind and sea. Ship cruising speeds are 
also usually constant and are determined by economical considerations, 
with slight variations to permit favorable entry conditions at the 
destination. 

There are exceptions to this rule of consistency. Fishing vessels 
off the Grand Banks follow very erratic maneuvering patterns. Navy ships 
executing training maneuvers or on certain patrol missions do not follow 
consistent course patterns. Since these exceptions apply to only a 
small fraction of world ocean traffic, a sequence of data points, i.e., 
position reports, which indicates a consistent course and speed is 
probably an actual ship's track. 

There are certain geographical areas where the antithesis of con- 
sistency is true. For example, ships approaching the Straits of Malacca 
from the north-east will, with a very high degree of certainty, change 
course radically to the north-west at the mouth of the strait. A highly 
sophisticated correlation program will be cognizant of the conditions 
and locations where these exceptions occur. However, even lacking this 
sophistication, much of the confusion resulting from circumnavigation 
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course changes will be resolvable by lending more credence to speed 
consistency than to course consistency in these areas. 

2.2 Application 

Each of the selection rules discussed have certain application 
limitations. One must conclude that an optimal correlation procedure 
would use them both to a certain degree, or in certain cases, but be 
able to distinguish when the selection rules break down. It is also the 
case that the longer the sequence of positions to be considered by the 
correlation logic, i.e., the smaller the interval between orbits, the 
more confidence one may place in the results. Certainly, there can be 
no less than a sequence of three points for consideration in order to 
initiate and extend tracks. 

The two correlation programs considered by the authors are described 
in detail in the next two sections. These programs employ both of the 
selection rules to some extent, but each emphasizes a different aspect. 
The first program, the deterministic correlation model, examines several 
sets of position reports at one time and makes maximum use of the selec- 
tion rule for initiating tracks, i.e., the elimination of consideration 
for correlation of two position reports from consecutive orbits which 
are separated by a distance greater than the distance attainable by a 
ship cruising at the arbitrarily stated maximum speed throughout the 
period between orbits. The second program, the probabilistic correlation 
model, emphasizes the selection rule for track extension, i.e., the 
greater reliance on consistent course and speed to select the most 
probable correlation. After considering each new data set, this program 
presents a list of all possible tracks together with a credibility factor 
which reflects the consistency of course and speed indicated by each. 

The probabilistic model will be described in detail in section four. 
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Section three contains a discussion of the logic, step-by-step proce- 
dure for correlation, and details the methods and results of the 
simulation used to evaluate the deterministic correlation model. 
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3. The Deterministic Correlation Model 

The deterministic model attempts to correlate reports of the 
position of ships at sea as obtained and relayed by surveillance satel- 
lites. Explicit use of probability theory is avoided. No reference is 
made to information such as normal cruising speed, or the existence of 
the more or less well defined shipping lanes between seaports. External- 
ly gathered intelligence information concerning shipping is not used by 
this model in its present level of development. The model is designed 
to permit track initiation and track extension by accumulating position 
reports from several orbits and correlating them concurrently. As it 
now stands, this model is a mid-ocean model. It has no logic for cor- 
relating position reports in proximity of seaports. 

To facilitate formulation of this model, certain assumptions have 
been made about the capabilities and characteristics of surveillance 
satellites and ocean shipping. 

These assumptions are: 

1) On every orbit the satellite borne sensor will report the position 
of all ships on the ocean that are within its range of detection. 

2) All positions reported will be positions of ocean shipping. 

Clouds, aircraft, land and objects on or over land will not be reported. 

3) Two positions reported at times differing by dt cannot corre- 
late if they are separated by a distance greater than the product of dt 

and S . where S is an arbitrarily stated maximum speed attainable by 

mx 9 mx 

ocean shipping. 

4) A ship reported on orbit i will also be reported on orbit i + 1, 
That is, the ship cannot enter port and thus avoid detection. This dis- 
regards the possibility of submarines being detected on orbit i and sub- 
merging before orbit i + 1. 
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3.1 The Correlation Logic 



Let P. be the set of n positions, p. , , 

i,n ’ *i,l ’ i,2 ’ ’ i,fi 



reported on orbit i. Let Q. . be the set containing all of the positions 

1 ) J 
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off set P. +1>m , which are within a distance d of p^ ^ , where d is the 
maximum distance a ship can travel in the time between orbits i and i + 1 



at orbit i will also be reported at orbit i + 1 and none will have travelled 
a distance greater than d, the maximum attaihable distance at speed . 

The step-by-step procedure for the concurrent correlation of ship 
position reports in the deterministic model iss 



2) Determine which, if any, of the Q. , have only one element, 

1 > J 

3) If any position which is the only element of a one element set 



4) Repeat steps 2 and 3 until no more changes can be made. 

5) For each j containing only one element a firm correlation 
exists between a position of orbit i and a position from orbit i + 1. 
Confirmation of this correlation is sought by projecting the track to the 
next (or prior) orbit. If the track is confirmed by the presence of a 
ship in the area of the projected position, the correlation is designated 
as a track. 



at speed , and p^ ^ is an element of P^ n 

of elements in the set Q. . . Then 

« J 
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Define n as the number 

q 




1) Construct the sets ^ for all p^ ^ 



is also included in any other set j 



remove it from those other sets. 
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6) Attempt to trace each designated track from the time of its 
designation back to orbit one and forward to the last orbit reported. 

7) Attempt to correlate tracks for the sets Q. containing more 
than one element. 

8) Attempt to correlate all remaining uncorrelated position reports. 

Step 6 above is accomplished by the following methodology: 

. a) Using the course and speed indicated for the designated 
track, determine the next projected ship position by projecting the track 
ahead (or back) to the time of the next (or previous) orbit. 

b) Define an area about this projected position and determine 
whether an undesignated position report from the next (or previous) 
orbit is located within this area. The size of the area defined can be 
varied as a function of the expected inaccuracies of reported positions. 

c) If only one ship position report is in the defined area, 
and if that position report is not simultaneously in the area defined 
about another projected position, it is assigned as an extension of the 
designated track. This position report is then removed from all other 
sets, Q. . , which may contain it. 

d) Repeat (a) through (c) until each designated track is 
extended as far as possible. 

e) Return to Step 2 until no further designations can be made 

on the basis of the deletion of elements from the sets Q. . . 

1 > J 

If more than one ship position report, or if no ship position 
report, is in the defined area, the projected track cannot be extended 
by the methodology for Step 6. 

At the completion of Step 6, the methodology for Step 7 is: 

a) For each Q. . consider the relation between p. . and 
1 1 J J 
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each element of j as a tentative track. 

b) Determine the projected positions for each of these tracks 
by projecting each tentative track to the next (or previous) orbit. 

c) Define areas about each of these projected positions. If 
only one of the areas contains a position report, this position report 
is considered as defining a designated track. 

d) For all tracks so designated, return to Step 6. 

Those Q- • which still contain more than one element, as a result 
of the areas constructed about the projected positions containing other 
than one ship, are analyzed to determine which ship is closest to the 
projected position. That position report which is closest is considered 
to be the most likely correlation and is designated. The track is then 
extended as in Step 6. 

In summation, the correlation logic for the deterministic model is 
based upon the previously stated assumptions and upon the following 
concepts: 

1) Any set Q. . which contains only one element determines a unique 
correlation with certainty of being correct. 

2) Tracks which are tentatively extended to a constructed area 
which is occupied by more than one ship, can be correlated on the im- 
plicit probability that the ship nearest to the center of the area is 
the most likely true correlation. 

3) Any track projecting to an area which contains no position re- 
ports does so as the result of a ship course or speed change. Based upon 
the implicit probability that a small course change or a small speed change 
is more probable than a larger course or speed change, the undesignated 
position report nearest to the center of the constructed area is considered 
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to be the most likely true correlation. 

Although some of the correlated positions may have been correlated 
erroneously due to changes in the course or speed of the ships, the 
assumption that each ship is reported on each orbit assures that event- 
ually all positions will be correlated. 

3.2 Applications of the Model 

A simple situation illustrating one of the basic principles of this 
logic is diagrammed in Figure 2. 




Illustration of One of the Basic Principles 
of the Deterministic Correlation Model Logic 



Figure 2 



Four position reports from orbit one and four position reports 



from orbit two are plotted and the Q. . have been constructed about 

J 

each P. .. The true tracks are drawn in solid lines with other poss- 
j 

ible tracks indicated by dashed lines. The position report p-, ^ can 
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To permit other 



be correlated only with members of the associated Q . 

1 > j 

correlations would violate the assumption concerning the maximum speed 
attainable. Due to the assumption that the position of each ship reported 
at orbit one will also be reported at orbit two, we are assured that each 
of the eight position reports plotted can be correlated. Therefore, the 
next step to be taken after the sets Q . have been constructed, is to 
determine which, if any, of the Q have only one element. From Figure 2 

1 9 J 

it can be observed that Q and Q have only one element. This deter - 

1,J 1,4 

mines a correlation between p and p and a correlation between p 

1,3 2,4 1,4 

and P 2>1 . 

It is clear that a specific ship cannot occupy more than one position 

at a given time. Therefore p n and p can be removed from set Q 

2,1 2,4 1,1 

and p_ . can be removed from set Q . Further examination now shows 

2,4 1*2 

that Q. p has only one element and a correlation can be established be- 

X,2 

tween p^ ^ and P 2 3 • After deleting p ? ^ from set a correlation 

is determined between p^ ^ and p^ ^ 

Although p^ , had four tentative tracks, a unique correlation of 

P, . and P . has been made. This simple, easily resolvable situation 

1,4 4,4 

is not likely to be encountered frequently, but it does serve as a basis 
and as a first step toward resolving more complicated situations. 

Figure 3 illustrates a situation requiring the appli ca ti on of all 
of the steps outlined in section 3.1 above. Four position reports from 
each of four orbits are plotted. Actual tracks are indicated by solid 
lines while tentative tracks and projections thereof are indicated by 
dashed lines. The circles indicate the areas constructed around pro- 
jected positions. The circles defining the sets Q. . are omitted for 

1 , J 

the sake of clarity. 
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The Q. are constituted as follows; 

J 

Q l,l = P 2,l ’ P 2,2 ’ P 2,4 

Q l,2 = P 2,l ’ P 2,2 ’ P 2,4 

Q l,3 = P 2,l ’ P 2,3 ’ P 2,4 

Q l,4 = P 2,l ’ P 2,2 ’ P 2,4 

The tentative tracks contained within Q have been omitted from the 

1,3 

figure for clarity. 

Since none of the Q contain single elements, the tracks are 

1, j 

initiated by projecting each of the tentative tracks and checking for 
position reports in the constructed areas. In this case only one of 
the tentative tracks projected from each of the p results in an 

1, J 

occupied area. Therefore all four tracks are initiated correctly. When 

track one and track three are projected from orbit three to orbit four, 

a conflict situation occurs. Both projections find an area occupied by 

two ships. This situation is resolved by comparing the distances of the 

tentative projections from the center of the area projected by the course 

and speed of the designated tracks, i.e., track one is projected to orbit 

four and finds p and p both in the area. However, p is nearest 

4,1 4,3 ’ 4,1 

to the center of the area and is the favored correlation. 

Tracks two and four, when projected from orbit three to orbit four, 

find only one ship in their respective projected areas, but it is the 

same ship, p. . A comparison of the distance from p to the center 
4 ,«£ 4 ,<£ 

of each of the projected areas shows track two to be favored and leaves 

track four with an empty area. The program must then look to the other 

members of Q, .. Since, due to the deletion of elements from Q q , as they 
->,4 3,4 

were assigned designations by other tracks p. . is the only remaining 

4,4 

element in Q, . it is correlated with track four. 

-*,4 
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Diagram of Track Initiation and Extension 
F igure 3 
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These two illustrations by no means point out all of the com- 
plications that may arise in track correlation. They are intended 
merely to demonstrate some of the fundamental points of the logic. 

3.3 Model Evaluation 

This correlation model was programmed for evaluation on the CDC 
1604 computer located at the U.S. Naval Postgraduate School computer 
facility. The program was written in the FORTRAN 63 computer language. 
The FORTRAN 63 program is listed in Appendix I. 

The program was designed to have the capability to receive a wide 
variety of inputs. Data representing specially designed shipping 
situations can be entered into the program. Alternatively, situations 
can be generated by the computer program by use of random number genera- 
tors. The program is written to correlate positions gathered from 
several orbits concurrently with none of the tracks designated prior to 
the initial orbit. It is considered that ability to resolve such a 
situation is more demanding than a situation where the program is merely 
asked to extend tracks previously designated and correlate them with new 
position reports. 

Parameters that can be varied to permit full analysis of the corre- 
lation program capabilities includes 

1) The number of ships to be tracked. 

2) The number of orbits to be processed concurrently. 

3) The time between orbits. For correlation of positions in a 
concentrated area, the speed of the satellite passage is so large rela- 
tive to the speed of the ships being tracked and the area concerned, 
that it does not invalidate the evaluation to consider that all of 

the positions reported during one orbit over the target area have been 
taken simultaneously. 
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4) The size of the area containing the initial positions of all 
ships. This parameter determines the shipping density. 

5) The distribution of initial courses. 

6) The distribution of initial speeds. 

7) The distribution of initial positions within the playing area. 

8) The distribution of the number of course and speed changes to 
be made. 

9) The distribution of the magnitudes of course and speed changes 
to be made* 

10) The size of the areas to be constructed about the projected 
positions. If it is required that the area size be incremented to look 
for ships that have changed course or speed, the program specifies that 
both dimensions of the area will be doubled. 

Provisions are also included that permit bypassing the routine 
which looks for the sq,t,sQ^ j which have only one element. When this 
provision is bypassed, the tracks are initiated by projecting all 
tentative tracks ahead to the next orbit and checking for area occupancy. 
Decisions are then made on the basis of minimizing the distance between 
the center of projected areas and the correlated positions. 

The first evaluation of the model involved the correlation of ten 
position reports received per orbit for five orbits. The time between 
orbits was chosen as one hour. The initial positions of the ten ships 
were uniformly distributed in a 200 mile square. Initial speeds were 
uniformly distributed in the interval 5-35 knots. Initial courses were 
uniformly distributed on the interval 0-360 degrees subject to the con- 
straint that no ship was assigned a course which would prohibit that ship 

from being included in the Q. , associated with a position of some other 

1 » J 

ship. No course or speed changes were made. 
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This situation was correlated once using the entire logic and 

once omitting the check for Q. , which contain only one element. All 

1 > J 

positions were correctly correlated, designated and tracked for all 

orbits. There was no detectable difference between the evaluation with 

and without the check for single element Q. . . 

1 > J 

The above situation was modified by changing the area containing 
the initial positions from a 200 mile square to a 100 mile square, and 
the evaluation was repeated with identical results. 

Using the same initial positions, initial courses and speeds, and 
the same area, the situation was modified to provide large course and 
speed changes at the time of orbit three for three ships. Additionally, 
one ship was provided with a new course and a new speed on every orbit. 
The model was again successful in correctly correlating, designating, 
and tracking all ships thru all orbits. 

The model was next evaluated using randomly generated data under 
the following constraints: 

1) Ten ships were to be tracked. 

2) Five orbits were to be processed concurrently. 

3) The initial positions of all ships was uniformly distributed in 
an area by use of a random number generator. The size of the area was 
varied as a parameter. 

4) Initial courses were generated from a uniform distribution on 
the interval 0-360 degrees. 

5) Initial speeds were generated from a normal distribution with a 
mean of 15 knots and a standard deviation of 5 knots, subject to the con- 
straint that the minimum speed was 5 knots and the maximum speed was 35 
knots. 
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6) For each of the last three time intervals the number of ships 
to change course and the number of ships to change speed were generated 
from a uniform distribution on the interval 2-8. 

7) The determination of which ship was to make each change was made 
by generation of a random integer on the interval 1-10. 

8) The magnitude of the course change was generated from a normal 
distribution with mean zero and with a standard deviation of 0.50 radians. 

9) The magnitude of the speed change was generated from a normal 
distribution with mean zero and with a standard deviation of five knots 
subject to the constraint that the resulting speed must lie on the inter- 
val 5-35 knots. 

10) Given a designated or tentative track at orbit i, the initial 
area to be searched for the extension of this track is centered on its 
projection from orbit i to orbit i + 1. The area extends - 2T miles 
radially and ± 0.09 radians in azimuth, where T is the time interval 
between orbit i and orbit i + 1. The radial dimension of the area is a 
function of T, while the angular dimension is a function of the distance 
from the area center to the track position at orbit i, i.e., a function 
of track speed and T. Figure 4 illustrates this area definition; 



X 

T r/tcfc H 
orfci-t L- 1 



—X 
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Definition of Area Constructed Around Projected Positions 

Figure 4 



2.1 



It is to be noted that the number of course and speed changes 
permitted for this test significantly exceed the number to be expected 
under any but the most peculiar situations to be encountered at sea. 
Further, the standard deviation stated for the magnitude of course and 
speed changes is considerably larger than that which would be expected 
under normal conditions at sea. These larger values are specified for 
this evaluation to provide a more severe test of the capabilities of 
this model. 

Figure 5 summarizes the test simulations, the simulation results, 
and indicates the parameters which were varied. A sample simulation 
and the results of the correlation are contained in Appendix II. While 
the availability of computer time prohibited conducting sufficient 
simulations to permit a rigorous statistical analysis of simulation 
results, the number of simulations conducted did indicate some general 
relationships between the various parameters and the results. A com- 
parison of columns five and six of Figure 5 indicates that a major por- 
tion of the incorrect track extensions involved the extension to the 
position reports obtained on orbit one. This is attributable to the 
much greater shipping congestion at that time due to the initial posi- 
tioning constraints. A major drop in program performance is indicated 
for the evaluations with an orbit interval of 8.0 hours and an initial 
area of 50 miles square or of 100 miles square. This drop is explained 
by noting that for each of the 80 simulations conducted under these con 
ditions, each of the Qj j constructed about the position reports from 
orbit one contained all of the positions reported on orbit two. A modi 
fication to preclude any course or speed change until three orbits had 
been completed was made for 20 simulations at 8.0 hour orbit intervals. 
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This required all ships to maintain a steady course and speed for 
16 hours. Under these conditions the effectiveness increased from U9% 
to 72% as indicated in Figure 5. In all of the simulated situations, 
many of the ships that were not correctly tracked through all orbits, 
were correctly tracked through four orbits. 
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Column 1 


Column 2 


Column 3 


Column 4 


Column 5 


Column 6 


SIZE OF SQ. 
CONTAINING 
INITIAL 
POSITIONS OF 
ALL SHIPS 


AVERAGE NO. 
OF COURSE 
AND SPEED 
CHANGES PER 
SIMULATION 


NUMBER 

OF 

SIMULATIONS 


ORBIT 

INTERVAL 

(HOURS) 


PERCENTAGE 
OF SHIPS 
CORRECTLY 
TRACKED 
THRU ALL 
ORBITS 


PERCENTAGE 
OF SHIPS 
CORRECTLY 
TRACKED 
THRU LAST 
FOUR ORBITS 


100.0 


26.9 


42 


1.0 


94.0 


96.2 


50.0 


21.8 


44 


1.0 


82.0 


90.5 


50.0 


16.3 


22 


2.0 


85.4 


93.7 


50.0 


23.8 


40 


8.0 


28.3 


86.5 


100.0 


23.0 


40 


8.0 


49.0 


90.7 


*100.0 


22.5 


20 


8.0 


72.0 


93.5 


500.0 


29.8 


20 


8.0 


81.5 


90.5 



* For these simulations, all course and speed changes were made for 
the last two orbits. 



Summary of Test Simulations for Deterministic Model 

Figure 5 
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3.4 Conclusions 

The following conclusions are permitted by the results of the model 
evaluations discussed in section 3.3 above: 

1) Under the stated assumptions, ships that maintain a steady 
course and speed will be correctly correlated, designated, and tracked 
in almost every case. 

2) For a shipping density of 10 ships per 50 mile square or less 
and for an orbit interval of 2 hours or less, the model performs track 
correlation with a high degree of reliability and precision, even when 
a large number of course and speed changes are permitted. 

3) For an orbit interval of 8.0 hours the model correlates tracks 
with a high degree of reliability and precision after the ships have 
departed from the congested area imposed by the initial problem condi- 
tions. 
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4. The Probabilistic Model 

The probabilistic model attempts to correlate all position reports 
obtained from the latest sweep of an area with all previously uncorre- 
lated position reports and all previously established tracks in order 
to determine all possible track extensions. A track reliability factor 
is then assigned to each extension to express the probability that this 
extension is the correct one. In this model track extension is primarily 
a function of the selection rule of velocity consistency as discussed in 
Section 2 while the track reliability factor includes those factors 
which must be considered in making a probability determination. What 
these factors are and how they enter the computation of the reliability 
factor form the logic of the probabilistic model. In the discussion to 

follow, the notation developed in Section 3 will also apply, i.e., P. 

i,n 

is the set of all position reports, p, Pi p i n re P or,ted from 

the i n sweep. 

In order to determine the reliability factor used in this model 
the possible results of correlating any one position report into a track 
extension are examined. For any point p. . from the set P, and a 

i j j i. 9 n 

current track A, the set of all possible events in extending track A to 

the set P. is? 
i.n 



E. 



,i p. . is the true extension of track A, or 

A o .1 



E, s p„ . is not the true extension of track A« 

1 , j 



The probability of the event E^ occurring, PCE^), is 



P(E 1 ) = 1 - P(i x ) 
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Upon further examination the event can be seen to include the 
following mutually exclusive events: 



E 2 : Pi j the "true extension of some other track X, 

while some other point p. , is the true extension 

i,k 

of track A. 



E^ : p^ j is the true extension of some other track X 

and track A has terminated since last sweep. 

E^j p^ j represents a ship which has put to sea since 
the last sweep, with no previous track, and some 
other point p. , is the true extension of track A. 

Ej.: p. . has no previous track and track A has terminated. 

? 1 , j 

The probability of E^ occurring is therefore 



p(i x ) = p(e 2 ) + p(e 3 ) + p(e 4 ) + p(e 5 ) 



and 



P(E X ) = 1 - ( P(E 2 ) + P(E 3 ) + P(E 4 ) + P(Ej) ) (4-1) 

It should be noted that inherent in equation (4-1) is the assumption 
that all ships are detected and their positions reported on each sweep 
by the surveillance satellite. 

The value of P(E^) obtained in equation (4-1) should be the value 
of the reliability factor used in this model. However, none of the 
probabilities in the right side of equation (4-1) can be computed directly. 
The value of P(E^) can at best be estimated. In the sub-sections to 
follow, equation (4-1) serves as a guide for obtaining a best estimate of 
P(E 1 ). 
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4.1 Probability Based on Velocity Consistency 



This section sets forth the methodology for determining a probabil- 
ity of correlation based on the selection rule of constant course and 
speed. Basic to this methodology is the assumption that course change 
and speed change are independent events. This assumption is based on 
the observations made in section 2.1 concerning circumstances governing 
course and speed changes. 

If the probability density function for speed change is continuous 
the probability that a ship will change speed exactly k knots is zero, 
regardless of the value of k. In order to obtain a positive probability, 
a small margin of error, A x, must be introduced. If f s (x) is the 
probability density function for speed change, the probability that a 
ship will change speed by an amount k i /\x iss 



However, this computation will still yield a very small probability 
even when the speed change, k, is zero, the most likely value. In 
order to determine a more meaningful value, therefore, the probability 



the probability of a correlation which shows the most likely speed 




k+/\x 



k>A x 



of a correlation p,,, showing a speed change of k ±A X as compared to 



change, 0 ~A X > is us ed« 

This probability of correlation is computed as the ratio of 



these probabilities, i.e., 
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p s = 



k +A 3 



f g (x)dx 



k -A- 



+A x 



fg (x)dx 



-A 3 



However, for constant A x in both numerator and denominator, this ratio 
can be assumed to be: 



P* = 



f s (k) 

f 9 (o) 



(4.2) 



The probability of the correlation based on course change 

p-. , is developed in precisely the same way and if f (x) is the probabil- 
u u 

ity density function for course change then 



p e = 



f e (k) 

~e (0) 



(4.3) 



where k is some representative value of course change. 

The probabilities, p g and p^ , represent two independent estimates 
of the probability for a valid track between position reports as a func- 
tion of consistent velocity. A best estimate of this probability is 
some linear combination of the two. In this model it is assumed that 

this best estimate for a given course and speed change is p _ , where 

s & 



Ps + P0 



S Q 
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The development of p s q thus far applies only for those tracks 
that previously existed in the model for which a course and speed 
were known. New tracks may be initiated in the model by connecting 
position reports from the previous sweep of the area which cannot be 
correlated with new position reports from the current sweep. In such 
a case there is no track history from which to determine a value of 
course and speed change. The discussion in section 2.1 implies that 
a truncated normal probability density function exists for ships' 
cruising speeds. For a correlation between a new position report and a 
previously uncorrelated report, the probability, p r , is used in this 
model instead of p g ^ , where 

f c oo 

P ° = f C 

and f c (x) is assumed to be the truncated normal distribution with 

l_L c = 14 knots and Cf Q = 2.5 knots and where k is the speed indicated 

by the time and distance between the position reports. This is to say 

that in the track extension phase, i.e., the phase for extending 

established tracks, p g ^ is assigned to all possible extensions in an 

effort to determine the most likely, while in the track initiation phase, 

p is assigned to all possible correlations between new position reports 
c 

and previously uncorrelated positions. 

For the evaluation of this model the density functions, f s (x) and 
f^ ( x )» were assumed to be truncated normal with the parameter values 
Ms = °* <Js = 2 *5 for s in knots and /jLq = 0, (Jq = 24 for 0 in degrees. 

No statistical justification is offered for the use of these distributions. 
They were selected because they seemed to be reasonable in the experience 
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of the authors. The methodology, however, would be the same re- 
gardless of the distributions used. It may be that experimentation 
or thorough statistical analysis will yield distributions which provide 
more accurate likelihood. These distributions, however, appear to 
adequately represent the conditions of velocity consistency for ocean 
going ships. 



are : 



Using these truncated normal distributions the desired probabilities 

— — e'^l (f ) 2 ,/,jV 

p 5 = *M!L = e -^l§ 



/ 



t(%j2TT 






(4.2) 



and in a similar fashion 




where Q = the amount of course change, 
k = the amount of speed change. 

= the speed indicated between position reports. 



(4.3) 

(4.4) 



t = the truncation factor. 

4.1.1 Probability Based on Mutually Exclusive Events 

The factor V s q expresses a likelihood that a particular track 
extension is a valid one based on the consistency of course and speed. 
In the probability sense, the set of all p s ^for the extension of some 
track A cannot be considered as a mass function to determine which 
extension of A has the highest probability of being valid. This is the 
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case because, in general, the set of all p g ^ for track A will not 
sum to unity. However this situation may be overcome by assuming that 



all possible extensions of track A from R to P i+ ^ n and the associated 
p g 0 form a set of mutually exclusive, exhaustive events. Thus a likeli- 
hood factor that represents the probability of a track extension being 
valid in view of the k possible extensions can be computed as 

( Ps0 )j 



L ( A Pi+l,j) 



£1 <r »e ) i 



where 



L(Ap^ + ^ j) is the probability of p ^ +1 j being the valid 

extension of track A with respect to consistent 
velocity and all other possible extensions of 
track A. 

is computed as in section 4 . 1 . 
is the total number of possible extensions of 
track A. 



( V 



It should be noted that the relative magnitude of the p _ have 

s U 

been preserved in the construction of L(Ap^ + ^ ^). For examples if 

k = 1, L(Ap. ) = 1.0 and if k = 3 and these three possible track 
i +1 f J 

extensions are equally likely then L(Ap^^ j) = *333 f° r a ll three 

regardless of the value of p . 

s & 

This same development applies to the computation of the probability 
of a newly initiated track being a valid one based on the speed indi- 
cated and weighted by the number of other possible new position reports 
which could be correlated with the previously uncorrelated position. 

The only difference is the use of p c instead of p s ^ wherever the latter 



appears, 



32 



4.1.2 Probability Based on Cross Mutual Exclusion 

The factor L(Ap^ + ^ j) is developed by a consideration of all 
position reports j ^ roin se ^ ^ + i m which fall within the dis- 

tance of the last known position of track A. It is also possible 
that the position report j is considered as a possible extension of 

more than one track, say tracks A, B, and C. If this is the case then 
there is an associated p _ as computed in section 4.1 for each of these 
possible tracks. The arguments presented in section 4.2 for the develop- 
ment of L(Ap^ + -^ j) will again hold and a new probability factor L(p^ + -^ ^A) 
can be defined which reflects the velocity consistency of P^ + ^ j with 

respect to track A normalized over the other tracks for which p. , . could 

i+i , J 

be the extension, i.e., 






(p *e 



l 

L (p 

i=l 



sQ ^i 



where 



L(pi + i jA) is the probability of p^ + ^ j being the valid 

extension of track A based on consistent velocity 
and all other possible tracks for which P^ + ^ j 
could be the valid extension. 



l 



is computed as in section 4.1, and 

is the total number of track extensions containing 

p i+i,r 



At this point the difference between L(Ap^ + ^ j) and L(p i+1 jA) 

should be noted. L(Ap„ . ) is the probability that position report 

»"-L A J 

p ^ , is the valid extension of track A with respect to all other 

i+i, j 
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position reports that may be valid extensions of track A. jA) 

however is the probability that position report p. , • is the valid 
extension of track A with respect to all other tracks that may contain 

p i+i,r 

The following computations can now be made j 



L ( Ap i+l,j) 



the probability that some position report, 
say p^ + ]_ jj, other than position report p^ + ^ j 
is the valid extension of track A. 

1 - L(Ap i+1>j ), 



and 



L( Pi + i,j A> 



the probability that position report p^ + j j 
is the real extension of some track other 



than track A 
= 1 - L(p i+1 jA). 

Referring to equation (4-1), P(E ? ) can be computed as; 



P(E 2 ) = L(*P i+1> j) X L(p l+M A) 

4.2 The Probability of Track Termination and Initiation 

Neglecting the last three terms of equation (4-1), a good estimate 

of P(E^) based on velocity consistency and mutual exclusion with respect 

to both track A and p. _ . is given by 

It J., J 

P(E 1 ) = 1 - P(E 2 ) 

This approximation of P(E^) does not include the possibility that track 
A has terminated or the possibility that p^ + ^ j represents a ship which 
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has just put to sea. Therefore, the possibility that tracks may be 
initiated or terminated between sweeps must be admitted for every 
track at each track extension phase. 

Probability of track initiation . If A is the mean rate of 
leaving port for ships all over the world then l/j^ is the mean inter- 
val between ships leaving port anywhere in the world. The average 
number of ships leaving port in a time interval t is therefore t/^ . 

If N is the average number of ships at sea the probability that any 
particular one of them left port in the time interval t is given by 



L ^ P i+l,j) 



VA 

N 



t 

T\ 



where, if t is the time interval between sweeps, then L(-p, , .) is 
the probability that the ship represented by the position report 
p i+l j pu ^ sea b e ^ ween sweeps i and i+1. The justification for 
computing -j) in the manner shown above is the observation 

that ships put to sea in completely random fashion. This assumption 
is false at any particular port because the rate of departure depends 
on conditions of weather, tide, and daylight or darkness. However, 
over the world as a whole the effect of these factors at one port is 
assumed to be cancelled by opposite conditions at another. There are 
not an infinite number of ships as would be implied by a completely 
random distribution; however, there is a relatively large number of 
ships in existence compared to the number of ships at sea, and the 
turn around time as compared to voyage duration is relatively short. 
The effect is similar to an infinite supply of customers in a queue. 

“ In the evaluation of this model the parameter values of A and N 
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were chosen as .10 and 5000 respectively. 

Probability of track termination . In order to determine the 
probability of a track terminating between sweeps the authors assumed 
a truncated normal distribution to represent the distribution of 
voyage lengths of ships at sea. Using this distribution the probability 
that a track will terminate between sweeps i and i+1, L(-A) can be com- 
puted as 



where d is the total distance for track A from the point of initia- 
tion up to its position on sweep i+1. 

There may be some justification in claiming that the probability 
of termination should be conditioned on the fact that the track was 
extended on a previous sweep, i.e., that the probability should be the 
probability that track A has terminated before sweep i+1 given that the 
track A had not terminated before sweep i. This conditional probabil- 
ity was not used in this model because the premise that the track had 
not terminated before sweep i is not a known fact. In the framework of 
the probabilistic model the track's existence on the last sweep is only 
a probability. 

Based purely on the experience of the authors, the parameter values 
of mean voyage length fj. ^ = 3000 nautical miles and Cf^ - 1000 nautical 
miles were chosen to test this model. More accurate parameters could 
be derived from a careful analysis of shipping records. 

It should be noted that if this model were to be expanded to include 
a computation of the distance of a position from known sea ports then 



A 



2 




0 



dx 
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both L(-A) and L(p 1+ ^ j) would be zero if the distance from port 



were found to be greater than S 



ffiX 



The values of P(E 3 ), P(E ), and P(E^) in equation (4-1) can now 
be found by using the values of L(-A) and L(vtp i+1 j), i.e., 






P(E 3 ) = L (Pi+l s jA) x L(-A) 
p ( E 4 ) = x L(-P i+1; j) 

P(E 5 ) = L ^“Pi+i s j) x L (“ A ) 



4.3 The Track Reliability 

Utilizing the individual probabilities calculated in the preceding 
sub-sections, the complete formula used in estimating the probability 
that a particular correlation is a 'valid one, P(E, ), is 

P(E 1 ) = 1 - ( l(Ap i+1 j) X L{p 1+1 jA) + L (p i+ i jA) x L(-A) 

+ i? ' A PuI~i7 X LC-P i+1 j) + L(-P 1+1J ) * U-A) 

This estimate is the track reliability factor used in this model to 
determine both the most likely new tracks and the most likely extensions 
for established tracks. In the next section a detailed sample problem 
will be used to demonstrate how track reliability is employed. 

4.4 Application of the Model 

The probabilistic model was programed in FORTRAN 63 to be run on 
the CDC 1604 computer. Summary flow charts and a program listing are 
included in Appendices III and IV. The logic of the program can best 
be shown using the set of diagrams in figures 6 through 11. 
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For illustration purposes assume that in the initial orbit of 

the satellite, the set of position reports P, _ is obtained. Since 

1,3 

there are no previous points with which this set can be correlated, 
they are retained in the model as uncorrelated positions. After the 
second sweep of the area a new set of positions P _ is recorded. It 
is determined that each of the points in the set P 2 3 is within the 
radius S^ of each of the previously uncorrelated points in set Pj 3* 
Figure 6 shows the nine possible tracks which must be considered. 



1,1 



''7' P 2,l 



A- 



1,2 x 






-Y P 2,2 






v' 

Y 






\ \ 

\ \ 

\ \ 
\ v 



1,3 



2,3 



Initial Tentative Tracks 
Figure 6 



A value of P(E^) is computed for each of the nine possible tracks 
and all nine are retained until the next sweep. 

At the third sweep a new set of position reports P^ ^ are received 
and the selection rule of consistent velocity can be applied. As before^ 
all positions from the set P^ ^ which fall within the radius of the 
positions in P~ ^ are considered possible extensions of the nine tracks 

2, j 
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established in sweep 2, and P^) is computed for each. In this track 
extension phase the probability p g ^ can be used instead of p c which 
was used in the previous phase. 




p i ( 3 f 2 ; 3 - 





Most Likely Tracks 
Figure 7 



In Figure 6 the most probable tracks are indicated by the solid 

lines. Utilizing the selection rule of consistent velocity to extend 

tracks the correlation (p-^ 3 P2 3 P3 3) has a higher track reliability 

factor than any other combination involving p, Hence (p, _ p 1 ) 

1,3 2,1 

and (pj - ? ) are discarded. By the same reasoning (p^ ^ p, ^) f 

^ P 1 1 p 2 3^ ( p ! 2 p 2 2^ and ^ P 1 2 P 2 3^ ar ® elijninated « Only three 
of the tracks produced by the position reports from sweeps 1 and 2 

| 

remain. 

The possibility that (p^ 1 P 2 2 ^ re P resen ts an initial track with 

a course and speed change to produce the track (p^ ^p 0 ^P-j 3) must also 

be considered when extending tracks from P.. _ to P_ The new possible 

1,1 3,3 
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track extensions from ^ to ^ to due to "near" constant 

velocity alone are illustrated in figure 8 by solid lines while those 
tracks possible due to a speed and course change are illustrated by a 
combination of solid and dashed lines. 




Tentative Extensions of Most Likely Tracks 
Figure 8 

For the tracks in figure 8 the solid lines receive a higher 
reliability factor than do the solid ~ dashed combinations. Figure 9 
illustrates the next phase of extension where the possible tracks are 
correlated with the elements from P. _. From the new set P the 
probabilistic model has assigned only two tentative track extensions 
which will be used in the next phase of track extension. 
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Firm Track Established 
Figure 9 

In the foregoing example tracks were initiated by the correlation 
of position reports from the current set P i+ ^ n with uncorrelated 
positions from the set P^ m . This type of track initiation will occur 
in the first few orbits of the surveillance satellite and may possibly 
occur at any time later. A much more likely situation in later orbits 
is that in which the position from the new set P i+1, n which represents 
a ship which has just put to sea will be considered as a possible ex- 
tension for one or more established tracks. It would not then be carried 
over to the next sweep as an uncorrelated position report. 



41 




Track Split 
Figure 10 



In figure 10 two tracks are shown, each with three possible ex- 
tensions. One of the points from the set P i+1, n must be a position 
report for a ship which has just put to sea. The decision as to 
which of the positions represents the new track is postponed until 
the next sweep. Figure 11 includes the position reports from the set 
?i+2 V computation of P(E. ) for each of the possible extensions, 

track 1.1 extended to p^ + ^ l determined to be the most likely 
extension of track 1, and track 2.3 extended to p^ + ^ j mos ^ likely 

extension of 2. This leaves the correlation - P ? +2 2^ aS 

new track. 
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New Track Initiated by Track Splitting 
Figure 11 

The track correlation problem for the probabilistic model as 
illustrated in figures 6 through 11 is, of course, for the particular 
sets, n# However consideration of these figures should give a 
basic understanding of the model. 

The correlation process is a matter of examining three consecutive 
sets of position reports to select a most likely combination of posi- 
tions. All but the most likely combinations involving the earliest 
position reports are eliminated. However, in going from one track ex- 
tension phase to another all combinations involving the last two sets 
of position reports are retained in the model. It should be noted that 
the elimination of the least likely backward extension is the equivalent 
of establishing the track's identity with respect to tracks previously 
established. If a track cannot be identified with a previous track, it 
must be assigned a new designation and considered to be a ship which has 
just come under surveillance. It should also be noted that retention of 
only the most recent position report together with a velocity vector is 
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the equivalent of retaining the last two position reports and a pre- 
dicted position for the next sweep. Thus, an identification symbol, 
a last known position, and a velocity vector provide the necessary 
and sufficient information for tracking in this model. 

4.5 Tests and Results 

In order to test this model a CDC FORTRAN 63 program labeled 
TRACK IV was constructed to conform to the logic of the model. This 
program and associated flow charts are included in Appendices III and 
IV. Another FORTRAN 63 program was developed in order to supply initial 
test data to TRACK IV. This input program is called READOUT and is in- 
cluded in Appendix III. 

The only formal tests imposed on the model involved variations on 
the situation shown in figure 12. While the results of these tests 
provide a basis for qualitative observations there was insufficient 
time to impose the variety and quantity of problems required for a 
quantitative analysis. 

The situation illustrated by figure 12 involved three tracks con- 
verging to an area and changing course when in close proximity. At 
about the time of course change for the three original tracks two new 
tracks were introduced. This problem was posed to TRACK IV on ten 
different runs with the interval between sweeps varying from one to 
twelve hours. Tabulated results of these runs are shown in table 13. 

'•Time to Resolve" is measured from time 64 when the problem became 
critical until the time when all tracks had been firmly identified and 
all course changes had been properly recorded. 

The column entitled "Tracks Correctly Extended All Orbits" applies 
to tracks properly initiated and correctly tracked without loss of 
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Tracks Correctly Changed Improper Time to Resolve False 
Interval Time To Resolve Extended All Orbits Identity Extensions New Tracks Tracks 
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Table 13 



identity throughout the whole of the problem. Improper tentative 
extensions may appear on any particular sweep with a higher track re- 
liability than the true extension. However, if on the next sweep the 
true extension is assigned the highest track reliability the track is 
considered to be correctly extended. It should be noted that the most 
common error which precluded a track from being listed as correctly 
extended was that at some point the track was redesignated. In such 
cases the track was never really lost for more than one orbit. Cases 
of redesignated tracks are tabulated under the next column, "Changed 
Identity". 

An entry is made under "Changed Identity" each time such error 
occurred. In the entries for the two hour interval one track was re- 
designated three times and another track was redesignated once. All 
of the errors which occur in the three hour interval concerned only 
one track. 

The column entitled "Improper Extensions" refers to cases in which 
a proper track jumped to assume the identity of some other valid track 
while the original track was either lost or assumed a new identity. 

This is considered to be the most serious type of error. 

The column heading "Time to Resolve New Tracks" refers to the two 
tracks introduced midway in the problem. The time is measured from the 
first appearance of the positions representing them until TRACK IV holds 
them as identified tracks with the proper extension being assigned the 
highest reliability. 

The last column lists those occasions where track designations 
were assigned to a group of possible correlations none of which were 
correct. On all three of the occurrences listed, the false tracks were 



4 ?' 



dropped on the next sweep. 

Another type of error was discovered in the tests which is not 
tabulated. In the one and two hour sweep interval runs* a number of 
tracks appeared which were duplicates of other firm tracks. This is 
due to a shortcoming of the program, TRACK IV, and is not a consequence 
of the logic of the model. The shortcoming is easily rectified by the 
addition of a sequence in the program which searches for duplicate 
tracks and eliminates the least likely. 

4.6 Recommendations a nd Conclusions 

The tests conducted indicate that the probabilistic model can be 
developed into a very useful instrument for assisting in the tracking 
of ships on the high seas. Only two of the shortcomings tabulated 
are inherent in the program logic. These are improper track extensions 
and the originating of false tracks. The latter is not considered 
serious because in each case the false track was in the system for only 
one sweep. Improper track extension occurs when at a particular read- 
out the inter-relationship between position reports indicates that two 
ships are making improbable course and speed changes. This fault can 
never be completely eliminated from the model but implementing some of 
the recommendations to follow may greatly reduce the frequency of 
occurrence. 

As a result of the tests conducted the following modifications 
should be made to TRACK IV before further tests are undertaken. 

1) More sophisticated programing methods should be employed in 
constructing the correlation lists. In the present model a maximum 
of 29 correlations can be tabulated for each track. The number of 
possible correlations increases as the square of the number of tracks 



in an area of congestion. Thus in its present form TRACK IV reaches 
a saturation level when 5 tracks converge 0 „ 

2) An internal track strength indicator should be assigned to 
each track. This can be a simple tally of the number of times a track 
is successfully extended. The track strength indicator will be used 
in connection with the next two recommendations. 

3) In its present form TRACK IV computes course change and speed 
change by taking the difference between those indicated by the last two 
positions and the last course and speed. This is accurate only if the 
ship changed course at precisely the time of the last position report, 

a very unlikely event. The result is that any course change appears as 
a simultaneous course and speed change neither of which are correct. 

This reduces the track reliability to a level far below what it should 
be, and probably is the cause of most cases of changed identity and 
improper extensions. A routine should be added to TRACK IV so that 
when track strength indicates that the track has been successfully ex- 
tended a minimum of perhaps three times and a course change is indicated 
with a simultaneous reduction in speed, the speed will be assumed to 
have remained constant and the exact amomt of the course change computed 
by solving a triangle with the law of cosines. 

4) Identity changes occur because of a shortcoming in that part 
of TRACK IV which eliminates the least likely backward extension of 
tracks previously listed as tentative. It is conceivable that all 
possible extensions for a track appear with a higher reliability as 
possible extensions of other tracks. In this case all the possible 
track extensions carrying the proper identity are eliminated, and the 
proper extension of that track appears as a split for some other track. 
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It is then handled as a new track. This situation can be remedied by- 
inserting a change in the sequence which handles track splitting so 
that the newly initiated track is compared with all tracks which have 
been dropped to see if it is in fact a continuation of one of them. 

If so and if the strength indicator of the old track was greater than 
three then the old track's identity and history should be adopted for 
the new one. 

5) Track duplication occurs as a result of several procedures 
in TRACK IV. The best way to rectify the situation is to add a sequence 
at the point in the program when all new track histories have been 
prepared. This routine would check for duplicate firm tracks and retain 
only those with the highest track reliability or those with the highest 
strength indicator. 

With the implementation of these recommendations the probabilistic 
model will be capable of handling much larger and more complicated 
problems and will yield much more reliable performance. 

The most surprising result from the series of tests of TRACK IV 
is the effect of extending the sweep interval. If those items used as 
headings in table 13 are used as criteria for performance then the ten 
hour sweep interval appears to be the most efficient from most stand- 
points. However, the results of the runs for sweep intervals between 
two hours and eight hours indicate that the factor which affects 
performance is the particular relationship existing among position 
reports on successive sweeps rather than length of the sweep interval 
alone. The model will have to be tested with a large number of varying 
problems before any significant observations with respect to sweep 
interval can be made. 
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5. Recommendations and Conclusions 

The following recommendations and conclusions are based upon the 
analysis of the evaluation of the proposed models. 

5.1 Recommendations 

The validity of the correlation models appears to warrant their 
further expansion and investigation to determine: 

1) The effect of relaxation of the assumption that all ships on 
the ocean and within range of the sensor will be detected and reported. 
In particular, further study should attempt to ascertain the lowest 
blip-scan ratio that will permit an acceptable level of correct corre- 
lati ons. 

2) The optimum size of the areas constructed about the track 
projections. Statistical analysis based on many random simulations 
should permit this determination. 

3) The maximum shipping density (number of ships per unit area) 
that will permit an acceptable level of correct correlations under 
various constraints regarding the number of course and speed changes and 
orbit interval. 

4) The effect upon the maximum shipping density of the input of 
externally gathered information concerning the course and speed of a 
specific ship. 

5) Whether these models can accept data gathered and preprocessed 
in the manner of the Pacific Missile Range simulation and correlate the 
data correctly. 

6) The appropriate measures of effectiveness to be applied to this 
model to facilitate analysis of the meaning and effect of various para- 
meter changes. 
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7) The design and effectiveness of a new model combining the 
essential features of the models proposed in this thesis. 

5.2 Conclusion 

The results of the evaluation of the two models proposed in this 
thesis indicate that track correlation using satellite gathered infor- 
mation can be satisfactorily accomplished by these and by similar 
correlation models. It is anticipated that further refinement of these 
models would provide a track correlation logic with very high capability 
and would be of significant value in further studies of the feasibility 
of ocean surveillance by satellites. 
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APPENDIX I 



FORTRAN 63 PROGRAM- DETERMINI STIC MODEL 

This appendix contains the computer program which was written to 
permit computer simulation and evaluation of the deterministic model. 
The program is written in the Control Data Corporation FORTRAN 63 
computer language and the model was evaluated by computer simulations 
on the CDC 1604 computer. An example of the problems generated by the 
program and a discussion of the input and output associated with the 
program are contained in Appendix II. 

Comment cards are used extensively in the program to explain the 
flow of the problem and to define the variables used. Detailed flow 
charts of the program logic are on file and may be obtained from the 
Chairman, Department of Operations Analysis, U.S. Naval Postgraduate 
School. 

The program uses the following FORTRAN 63 library functions: 

1) RANF 

2) SINF 

3) COSF 

4) ATANF 

5) SQRTF 

The program compilation time is 2 minutes and 53 seconds. Com- 
putation time ranged from 8-18 seconds per simulation. 
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ORBIT N TO POSITION J FROM ORBIT N+l. 

B ( I » J ♦ 3 ) CONTAINS THE COURSE REQUIRED TO GET FROM POSITION I FROM 
ORBIT N TO POSITION J FROM ORBIT N+l. 

B ( I ♦ J * 4 ) CONTAINS TRACK DESIGNATION FOR CORRELATED POSITIONS. 
SIMILAR FORMAT HOLDS FOR PHASES 2, 3» AND 4. 
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GENERATE INITIAL SPEED FROM NORM ( 1 5 , SPDDEV ) , SUCH THAT SSPD GE 5 AND LT 35 
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1*1=1* NPOS I T ) 

2001 FORMAT (/♦11X*1HT»I2*8X»F8«3»6X*F6.3»10X» 1 1 , 1 2X *F6 . 1 » 12X » F6 . 1 ) 



DO 1004 JA= 1 » NORB I T 
DO 1004 I A= 1 * NPOS I T 
A ( IA,JA»1)=JA*TIME 
A ( I A » JA » 2 ) = SX ( I A * J A ) 
1004 A( IA»JA»3)=SY( I A * JA ) 
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COMPUTE DISTANCE, COURSE AND SPEED FROM EACH POSITION OF ORBIT 
EACH POSITION OF ORBIT J+l FOR J= 1 »NORB I T- 1 . ENTER IN B(IB,JB, 
IF SPEED IS LESS THAN OR EQUAL TO MAXSPD. 
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PROJECT TO NEXT PHASE TO DETERMINE NUMBER OF AREAS OCCUPIED 
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IF ONLY ONE AREA FOUND OCCUPIED BY COLUMN CHECK, THEN CHECK THE 
ROW CONTAINING THAT ELEMENT TO SEE IF ANY OTHER ELEMENT COULD 
PROJECT TO SAME AREA. IF SO, CANNOT RESOLVE, GO TO NEXT COLUMN. 
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REMOVE ENTRIES FROM ROWS AND COLUMNS OF C AND D FOR THOSE DESIGNATED 
CORRELATIONS. 
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503 I B=N $ IF(B( IB, JB,L+4) .GT. 0.0) 500,504 

504 NUN I QC = 1 $ CALL PRO JECT ( I B , JB , MP ) $ NUNIQC= 
500 CONTINUE 

RETURN 

END 



SUBROUTINE UNIOROW(MP) 
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SUBROUTINE PRO JECT ( I B » JB * MPP ) 
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IF ONLY ONE POSITION IN AREA, THEN CHECK TO SEE IF ANOTHER TRACK 
CAN PROJECT TO SAME POSITION. IF NOT, DESIGNATE. 
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612 I P =NB I NR S JP = NB I NC 
GO TO 613 



IF PHASE GREATER THAN ONE» PROJECT BACK TO EARLIER PHASES 









IT 







o 


xO 


o 


r-H 






xO 


o 


r*H 


r-H 


CM 






r-H 


r-H 


vO 


CM CM 


xO 






vO 


vO 


» 


nO xO 


* 






•k 


* 


On 


♦ * 


lT\ 






00 


CN 


r-H 


vt in 


CM 






r-H 


r*H 


NO 


CM CVJ 


xO 






V0 


vO 


4— * 


xO vO 




r-H 






*"■* 


u 


N 


U 


1 




X 


X 


X 


X X 


X 


< 




1 — 


t— 


h- 


f— h- 


►— 


CL 




CL 


o 


O 


CL O 


o 


21 


o 


UJ 


t— 4 


t— * 


UJ •— 


»_4 


II 


il 


O 


3: 


> 


o y: 




< 


2 


CO CO 


CG 


CO CO 


CO 


CL 


UJ 


• 


• 


• 


• • 


• 


2 


_J 


UJ 


UJ 


UJ 


UJ UJ 


UJ 


** 


UJ 


_J 




o 


-J _J 


o 




2 


• 


• 


• 


• • 


• 


r-H 













*■— 


4— • 


«— k 






























4— 










II 
















— -* 


— 








































u 
















— 


<■— 































<— > 


^ — 








0 




■> r-H 










CM 


r-H 


r-H 




























CM 


cn 


cn 












II 










1 


1 


1 




























-f 


+ 


+ 








2 




o 










-J 


-J 




u 


vO 


















O 




r-H 


_J 


— J 


-J 






in 


O 




< 






xO 




»• 


#k 


•k 


t-4 


CO 








On 




vO 






r-H 




CM 


•k 




•k 






rH 


2 




_J 






rH 




U 


u 


u 


II 


vO 








cn 




CM 






CM 




xO 


0 


O 


O 






xO 






Ll 






xO 




1—4 


»-4 


♦—4 


u 










xO 




xO 






vO 




•k 


*— 4 


1-4 


1-4 






•k 


& 










•k 




•k 


•k 


•k 


2 


O 








•k 










•k 




cn 


#k 




* 






xO 


^k 










m 




CL 


CL 


CL 


♦—4 


CM 








rH 




0 






CM 




CM 


UJ 


UJ 


UJ 






CM 


U 




fcA 






cn 




“) 






CO 


xO 












►— 






CM 




x£) 


*-4 


>-4 


*-4 






xO 


2 










xO 




—• 


— 


W 


2 










xO 










xO 






W 


w 


w 






- — . 


t-4 




CO 










CO 


CO 


CD 




r-H 












0 












CD 


CO 


GO 






CM 


CD 




-> 






0 




1 


1 


1 


CL 


• 




m 




O 




0 






O 




0 


1 


1 


1 






• 


2 




II 






• 






— 




“) 


O 


m 


r-H 




• 










• 




• 


«— 


r—» 


<— . 






O 


•k 




m cl 






0 




CM 


CO 


CO 


11 


UJ 


rH 


xO 




O 




b9 






O 




0 


CM 


rH 


rH 






UJ 


CL 


r- 


r-H -) 






• 




+ 


+ 


+ 


OC 


• 


vO 


#k 




• 




U 




1— 


• 




• 


1 


1 


1 






• 


1-4 


CM 


xO 




h- 


I— 




-J 


_J 


—1 


2 


2 


•k 


00 


l— 


K— 




•—4 




»—• 


1— 


1— 


h- 


_J 


-J 


-J 






2 




xO 


* ^ 




t-4 


O 




•k 


#k 


•k 




»-4 


r- 


cn 


*— 4 


O 




<1 




LO 


KD 


1-4 


O 


•k 


*k 


* 






«— 4 


CL 


•k 




1 


if) 


• 




CL 


CL 


CL 


CO 


CD 


co 


xO 


if) 


• 




u 




O 


• 


CO 


• 


U 


u 


u 






CO 




m 


r-H CO < 


0 








-> 


-) 


2 


2 


vO 




O 






2 




CL 




O 


— - 


•— 4 










2 


•k 


rH 


xO 


CL 


CL 


cn 




•k 


•k 


r> 




W 


— 


r-H 


CL 


cn 




*— * 




2 


cn 


CL 


rH 


•k 


•k 


•k 






w 


< xO 


II 


2 


2 


1 


r-H 


CL 


CL 


CL 


b9 


Ll 


r-H 


• 


2 


1 




CO 




♦. 


1 


2 


+ 


O 


a 


Q 






Ll 


CL 




— CL 


* 


•> 


_J 


+ 


r — 4 


►-4 


M 




*—• 


• 


O 


•k 


-J 




2 




rH 


— J 


•k 


_J 


1-4 


*-4 


*— 4 








2 


<— • 


r-H M 




r-H 




2 


w 


w 


w 


r-H 




O 


UJ 


r-H 


•k 








II 


«k 


rH 


•k 


w 


w 


W. 


r-H 






w 


rH 


• 


11 


II 


u 


UJ CO 


CO 


CO 


+ 




UJ 


• 


II 


u 




& 




O 


u 


II 


Q 


CD 


CD 


CD 


+ 




& 


0 


• 


1- 


-J 


u 


K- 4 


-J 


w 


w 


w 


2 




• 


u 


u 


»— 4 










t— 4 


UJ 


»— 1 


w 


w 


w 


2 






*— 4 


O 


o 




1-4 


#k 


UJ 


u_ 


u. 


U_ 


1-4 


UJ 2 


0 




•k 


UJ CL 


u 




•k 


t-4 


#k 


ll 


U- 


Ll 


•-4 


UJ 


UJ 


f) UJ 


• CL 






CL 


2 tO if) 


if) 


CD 


3 UJ 


—4 




CL 




2 


O 


O 




UJ 


O) 


f) 


(■ OCDO 


Z> UJ 


• 


CL CL 


O 


xO “) 


II 


CO 


CO 


CO 


2 


2 


_J 


2 


O' 


“> 


2 


II 


1-4 


rH 




rH 


*-4 


□0 


CD 


CD 


2 2 2 0 


< 


a 2 : 


11 


r-H 


w 


2 < 


< 


< 


II 




UJ 


X CO 


4-- 


•— 4 


(X 


CO 


CM 


4 *— 


CM 


w 


< 


< 


< 


II 




♦—* 




CL 


2 : 11 


2 


NO CO 


UJ 






'k-' 


2 


h- 


2 


2 


xO CO 


1— 


2 


2 


xO CO 


xO 


CO 


W 


W 




2 


h- 


h- 


-J 


2 


w < 


♦— « 




w 


-J 


w 


— 


W 


*— • 


2 




— 




—4 


2 


♦— 4 


II 




w 




4— 


W 


W 




♦—* 


2 


2 


-J 


w 


U- CL 


CO 0 U- UJ 


Ll 


u_ 


Ll 


CD 


O u- 


Ll 


O Ll 


O CO 


u 


O u. 


0 u. 


ll 


Ll 


U.C0OO 


< 


Ll 


— 2 


2 


0 


t— 4 


2 


— 


*~4 


*— 4 


2 


u 




•“4 


O 


*— 4 


u 


2 


*— 4 


Q 


— 


0 


I—* 


»— 4 


*— 4 


K-4 


2 U U 


u 


1—4 


m 


00 






m 


r- 


CO 


0 


On 


xO 


xO 


r- 


00 




On 


rH 


O 






CM 




cn 




0 


m 


r-H O 


xO 




O r-H CM 






CO 


r-H 


r-H 


On 


r-H 


r-H 


co 


co 


CO 




cn 




CM 






CM 




CM 


CM 




CM 


CM 


r-H 


CM 




xOxOxO 






xO 


xO xO 


r-H 


vO 


vO 


xO 


vO 


xO 




xO xO 


xO 






x£> 




xO xO 


CM 


xO xO CM 


xO 


















vO 
































xO 






xO 







U U ( 



76 



627 IP=NBINRSJP=NBINC$NELEM=0$ GO TO 628 
615 NUN I QC = 0 $ RETURN . 

END 
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757 TR ACK ( N »MPA-1 ) = I C $ B ( JB , I C , L ) =N$ B { I B , JB , L+4 ) 
TRACK(N»MPA)=JBSTRACK(N»MPA+1)=IBSG0 TO 730 
707 RETURN 
END 
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APPENDIX II 



INPUT AND OUTPUT - DETERMINISTIC MODEL 

The variable names used by this program are defined in the program 
listing contained in Appendix I. No provision is made for reading the 
values of parameters into the program. Changes in parameter values are 
made by inserting punched cards into the program deck prior to compila- 
tion. The parameters which determine the formulation of the problem 
are : 

1) NPOSIT 

2) NORBIT 

3) SI2ESQ 

4) CUSDEV 

5 ) SPDDEV 

6) MAXSPD 

7) MINSPD 

8) MX 

9) NUNIF 

10) NPROB 

11) BDEPTHO 

12) BWIDTHO 

13 ) BWIDTHCO 

14) TIME 

The DIMENSION statement for all arrays must be changed if the value of 
NPOSIT exceecfe 10 or if the value of NORBIT exceeds 5. The values of the 
other parameters may be changed freely. 
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The sample output which follows illustrates several points: 

1) Governing parameters are printed at the beginning of the out- 
put. The number of course changes and the number of speed changes to 
be made for orbit J, as determined by random numbers, are indicated. A 
course change or speed change indicated for orbit J determines the course 
or speed for the period between orbit J-l and orbit J. 

2) The column labeled TRACK gives the position of each ship at the 
time of the orbit number listed in column TIME. The course and speed 
listed is the course and speed for the ship indicated for the period 
prior to the time of the position report. 

3) The course, speed and position for each ship at the initial 
orbit is generated by use of random numbers. The program then determines 
how many changes of course and speed will be made, which ship will make 
them, and the sign and magnitude of the change. The later courses, speeds, 
and positions are then computed and used as a basis for correlation. 

4) NCYCLE indicates the number of times the dimensions of the 
constructed areas were increased before the final correlation was made. 
Final dimensions are indicated. 

5) The print-out of tracks indicates which ship was designated as 
a specified track for each orbit. Any track which contains the same 
ship number for all orbits is considered to be perfectly correlated. In 
this example Track 2 and Track 14 were erroneously projected back to 
orbit 1. The cause is not apparent without plotting all tracks. 

6) Note that a "mistake" in correlation as indicated in this 
example does not necessarily imply faulty correlation logic. Many cases 
of this type result from crossing tracks and erratic ship maneuvering. 

The "correct" logic may pick the "wrong" ship. 
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7) The tracks containing all zeroes result from duplication of 
designations. For example, ship number seven may have been designated 
as Track 3 for orbits 3, 4, and 5. It may have been later designated 
as Track 6 for orbits 1, 2, and 3. The duplication is resolved prior 
to printing the results of the correlation. 
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APPENDIX III 



FORTRAN 63 PROGRAM - PROBABILISTIC MODEL 

The input to TRACK IV is a magnetic tape prepared by some program 
such as READOUT. This input tape may contain simulated position re- 
ports for one or more readouts. The data from one readout can contain 
from one to 100 position reports and must be written with a FORTRAN 
'WRITE TAPE i,L' statement where L is an array dimensioned as (100,3). 
Each position report occupies three words of the array L as follows: 
L(i,l) = Time of report 
L(i,2) = Latitude of the position 
L(i,3) = Longitude of the position 
The array L must be TYPE REAL, that is all data must be in floating 
point. The data must be in the units of measurement described in the 
next paragraph. A problem termination is indicated on the input tape 
by an END OF FILE. No END OF FILE should appear between readouts. 

The method and units of measurement for the quantities processed 
by TRACK IV are as follows: 

1) Time - Time is measured in hours from the beginning of the 
problem to the end. Time need not begin with 0.0 for the first report. 
In fact, if the program READOUT is used a time less than 2.0 will not 
be processed. Time is a floating point quantity. 

2) Latitude - latitude is measured in radians from 0.0 radians 
at the north pole to 77 " radians at the south pole. It is a floating 
point quantity. 

3) Longitude - Longitude is measured in radians from 0.0 radians 
at the Greenwich meridian eastward through 2 TT radians. It is a float- 
ing point quantity. 

4) Course - Course is measured in radians from 0.0 at north to 
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the right through 2 7T radians. It is a floating point value. 

5) Speed - Speed is measured in radians per hour. It may be 
converted to knots by multiplying the value by 3436.6. Speed is a 
floating point quantity. 

6) Distance - Distance is measured in radians and may be con- 
verted to nautical miles by the same conversion factor specified for 
converting speed to knots. It is a floating point quantity. 

7) Probability - Probability is handled in two forms. In most 
computations it is a floating point value between zero and one. It 
is stored in lists by multiplying the value by 1000 and converting 
to integer form. This permits the quantity to occupy 12 bits in 
memory. Some computations in TRACK IV handle probabilities in their 
integer form. 

TRACK IV utilizes a permanent record tape on which to store track 
histories between readouts. The tape also contains the next available 
track number and a list of any uncorrelated position reports. The 
format of these lists is as described below for the arrays TRACK and 
OLD DATA. This tape is set to an initial blank form by TRACK IV before 
processing the first input from the input data tape. If it is desired 
to continue a problem on different runs, that part of the program which 
initializes the permanent memory tape must be removed before the second 
run. 

Two tape units are used by TRACK IV as temporary storage. These 
are identified as units two and three in the program. The input tape 
is designated unit one and the permanent memory tape as unit four. 

Data design for TRACK IV is as follows: 

1) OLD DATA and NEW DATA are arrays dimensioned (500,3) and 
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(100^3) respectively. Both contain position reports in the form: 

NEW DATA (1,1) = time of report 
NEW DATA (1,2) = latitude 
NEW DATA (1,3) = longitude 

2) TR LIST and TR VEL are equivalent names for the same array 
which is dimensioned as (31,4). Data is addressed as TR LIST if it 

, is an integer value and as TR VEL if it is a floating point quantity. 

a) TR LIST(1,1) bits 0-5 contain the tentative part of 
the track number. Bits 6-23 contain the basic track number. Bits 
24 - 47 are unused. 

b) TR LIST(1,2) bits 0-11 contain the probability of 
termination computed on the last sweep. Bits 12 - 23 contain the 
track reliability assigned to the track on the last sweep. 

c) TR VEL(1,3) and TR VEL(1,4) contain the time of initial 
track contact and the time of last extension respectively. 

d) TR VEL(2,l) and TR VEL(2,2) contain the latitude and 
longitude of the last known position respectively. 

e) TR LIST(2,3) is set to 777 octal if the list concerns a 
previously uncorrelated position report instead of a previous track. 
This flag indicates that there is no history for the track. 

f) TR VEL(2,3) and TR VEL(2,4) contain the last known 
course and speed respectively if the list concerns a previously known 
track. 

3) The remainder of the TR VEl/TR LIST array contains possible 
correlation in the form; 

a) TR UST(I,1) = new data list index corresponding to the 
location of the position report in NEW DATA. 
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b) TR LIST(I, 2 ) p g 0 or p c for the correlation listed 
until the list is normalized? L(Ap.j^ ,) thereafter. 

c) TR VEL(I,3) and TR VEL(l,4) contain course and speed 
respectively as indicated if this correlation is a valid one. 

4) PL LIST is dimensioned (100,29) and is TYPE INTEGER? that is, 
data in the list is integer rather than floating point. In this array 
all correlations which involve a specific new position report are 
tabulated. Each entry requires one computer word. Bits 0-11 contain 

the p or p^a as appropriate until the list is normalized? L(p, . .A) 

C 3 (7 

thereafter. Bits 12 •= 17 contain the track tentative number and bits 
18 - 32 the basic track number for the track with which the position 
report may be correlated. Bits 33 - 47 are unused. An entry addressed 
as PL LIST(I,N) refers to the N ' possible correlation for the I entry 
in NEW DATA. 

5) TRACK NR and TRACK are equivalent names for the same list. 

The list is eight computer words long. TRACK contains the permanent 
record form of the track history as it is stored on the permanent 
memory tape. If a data word in the list is an integer it is addressed 
as TRACK NR? if it is a floating point value it is addressed as TRACK. 

a) TRACK NR(l) contains the track number and tentative 
number in bits 0—20. Bits 21 — 47 are not used. 

b) TRACK NR(2) bits C - 11 contain the probability of termin- 
ation computed at the time of the last extension. Bits 12 - 23 contain 



the track reliability. Bits 24 - 47 are not used. 

c) TRACK (3 ) to TRACK ( 8 ) contain, in the order listed, time 
of track initiation, time of last position, latitude of last known 
position, longitude of last known position, last known course, and last 
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known speed. 

Table 14 is a sample output from TRACK IV. The numeral to the 
left of the radix point in the track number column is the basic track 
number. Ideally it would remain unchanged from the time a track is 
initiated until the time it is terminated. The digits to the right 
of the radix point are branch identifiers. They are assigned at each 
readout to identify all possible correlations for the basic track. 

Arrows to the left are not part of the print out. They were 
inserted by the author to show which of the possible extensions is 
the true one. The sample shown was the output for the test described 
in figure 12 for the four hour interval run. 
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SEA SURVEILLANCE - PROCESSING TIME 68.000 
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Sample Output for TRACK IV 



-COOP, , MCARTHUR SURVEILLANCE , I / 1 /4 /0/4/S / 56/ 57 /£/ 2 = 56/ 3 = 5 7 
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INITIALIZE PERMANENT MEMORY TAPE 






o 

o 



LU 



II «— 



Z — 



LU 



3: X 
LU LU 
CX Z 



vO 

O 

O 

o 



r- 

o 

o 

o 



‘0106 (3 
8006 ‘600 




o z 

4—4 

^ H— 

o u 


LO 


• o 


o 


o cx 


o 


CL — 


o 


• «-H 

O • * 

liJ K ^ 





o — 



o 






















• 


< 








o 








• 


< 




• 


















1 — 








o 








o 


>— 




LD 


< 


<\J 








o 


r-H 




<C 








o 










< 




Z 


h- 


O 








o 


#N 


i — ♦ 


O 
















II 


p 




4—4 


< 


O 








r— H 













r-H 














l/) 


o 


O 






o 


«K 




*— i 


IS 
















— 


o 




u0 










• 


r-H 


< 


W 


LU 






l 


•— < 










_J 




LU 


3 








o 




1 


< 


Z 


















o 




U 


LU 










II 


C 


1 — 








*— i 


• 








t-H 




vf 


O 


z 


• — * 






II 




o 


<C 


II 








o 














QC 
















o 








ii 


LU 








< 


<r 


LU 


a 


*— » 


•> 


I—* 




LU 




3 




LU 


LU 






• 




CNJ 


CO 


1 — 




_J 




H 


Ll 






Z T- 


Ll' 


IS 




3 


r-H 


cn 










<c 






LU 




o 


o 




>—4 


o 


Z 


LU 


*— » 


Z 


o 




t—t 


o 


o 


o 


o 


LU 


LL 


u 




LU 


z 




h- 


o 




z 


h- 


4—1 


» — I 


X 




z 


z 


z 




I— 




z 


o 




*— * 


a 


U 


o 






U 


1- 




LU 


— > 


t—t 


►-HI 


4—4 


o 


»— « 


o 


LU 


< 




3 


O 


o 








o z 


II 


O 




3: 






u 


0- 


z 


Z 


LU 


Ll 


LU 


i— 


QC 


O 


Ll 


Ll 


cv O 




z 


Ll 


UJ 


LU 


LU 


o 


3 


LU 


Z 


X 


*— < 


QC 


t/) 


a 


O 


►— < 


*—4 


CL 


u 


*— 1 


» — * 


*— i 


QC 


X 


(X 








o 


































r— < 






u 


|H 




<\J 




vO 






00 


o 


r- 




o 




m 














o 




O 




o 






o 


r— < 


o 




o 




o 














o 




O 




o 






o 


o 


o 




o 




o 














o 




O 




o 






o 


o 


o 




o 




o 







u u u 



u u u 



94 



KEY 2 IS SET IF A STOP IS DESIRED JUST BEFORE WRITING FINAL TAPE 
KEY 3 IS SET IF A STOP IS DESIRED BEFORE PROCESSING NEXT NEW 
DATA 



00 

< 

Cd 



cc 

o 

U_ 

V 

u 

Uj 

X 

u 



o 

UJ 

o 



in 

o 



</) 

i— 

2 : 

O 

Q_ 

O 

UJ 



UJ 

Cd 

Cd 





0 


O 






O 


< 




0 


f\l 


O 




u 


h- 




rH 






2: 


2: 


< 




•s 




II 


h- 


Z> 


Q 


h~ 


r — 1 


1 — l 










00 






— * 


1— 


0 


O 


■ 4 


II 


II 


n 


X 


_J 


_J 


_J 








UJ 


0 


O 






~) 


*— 1 


2: 






1 — 1 






■ — ■ 




CL 


— » 


CL 


f — ( 


rH 


h- 


— - 


O 






rH 


r-H 


00 




u. 


■ — 


cd 


O 


O 


> — < 


«-r 






< 


rH 


rH 


_j 


0 


h- 


O 


UJ 








< 


00 


< 


-J 


O 


O 


i — 1 


UJ 


UJ 


UJ 


u 


O 


Q 


a. 


cd 


I— 


Cd 



00 

2: 

o 



LlJ 

Cd 

Cd 

O 

u 

UJ 
I 

CD 

00 

00 

O 

CL 

Cd 

o 

Ll_ 

V 

U 

UJ 

X 

u 



o 

o 

OvJ 



o 

UJ 



o 

o 





O 


1 — ( 






O 






It 


CO 


<r 


II 


O 














in 




h- 




» — f 


• — 1 




* 






— 






") 


r-H 




rH 


►— « 


II 


11 


e* 


11 




II 


w 








*— « 




•— H 


Cd 


•— < 


*> 


— 


«— ■ 






2 : 






1 — 






CM 




CO 


CO 


00 


<r 




O 




0 


0 


>— i 


0 


CO 


»— I 


u 


1 — ( 


rH 


-J 


rH 






< 










II 


O 


Cd 


0 


O 


CL 


O 




Q 


I— 


0 


O 


h- 


O 


*0 


»— < 


<\J 






CO 






O 


O 






0 






«— ( 


» — l 






rH 







< 

n 

o 

-j 

o 



m 

o 



u u u u u u u 



u V u 



u u u 



u u 



95 



LAT 1 = OLD DATA (11,2) 
LON 1 = OLD DATA (11,3) 
TIME 1 = OLD DATA ( 11,1) 
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COPY DATA IN PI LIST 
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120 TR LIST! J» I ) 



COMMENCE CHECKING TRACKS FOR CORRE LLAT I ONS 
200 INDEX 0=11-1 
2001 READ (A) TRACK 

IF (EOF, 4) 300, 201 
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PT IS IN VICINITY OF THI.S TRACK. SEE IF CORRELATION POSSIBLE 
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ALL POSSIBLE CORRELATIONS HAVE NOW BEEN LISTED. CHECK UPDATED 
TENTATIVE TRACKS FOR DUPLICATIONS. 
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ASSIGN 3017 TO NEXT 
3017 READ (2) TR LIST 
I F ( EOF , 2 ) A00 ♦ 301 

301 IF ( TR LIST (2.3) .EQ. 777B11000, 302 

302 IF (TR LIST (3.1) .EO. 0) 320, 330 



STEP 320 REPRODUCES A TRACK WHICH WAS NOT UPDATED THIS SWEEP 
EXACTLY AS IT WAS BEFORE. 

IF (PROCTIME - TR LIST(1,4) .GT. DEADTIME 13017,321 



V 







x 




h- 1 






LU 




Z 






UJ 




D 




» — i 


5 










tO 




UJ 




GO 






X 




o 


h- 




O 




(X 


to 




h— 




CL 


< 




to 














• 






o 






UJ 




Z 




o 


> 




< 




• 
















X 




GO 


< 




< 




o 


h- 








o 


Z 




_J 




o 


UJ 








o 


h- 




UJ 


CXJ 


r- 






1— 


cn 


co r- 


to 




X 


CO 


o o r- 


< 




a 




o o r- 








m 


o o 






o 


CO 


O 1 — * • 


X 




u 


.CO 


r- x O 


u 








r- ~ z 


►— 1 




•V 


— ■>> 


r** r-H < 


X 




UJ 


rH 


• 


5 




> 




CO 








CD 


• o — 


V 


GO 




O 


O X CM 


u 


c-* • — 1 


< 


(X 


^ x - 


< 


I s - co 


\— 


a 


<c ■— < 


x 


CO 


Z 




• * — 




• 


UJ 


• 





f— 


1— 


co 




1— 1 


(\J 


(O 












* 


i — 






1— 


— - 








to 


to 


to 


















< 




o 








o 


(NJ 




C\J 


<XJ 


_4 


*— ' 




i — i 


1 — 1 


f\J 


(XJ 


(XJ 


C\J 








< 


o 


>- 






• 








I — 1 


_J 




X 


W 


— 


— * 




— * 








X 


• 


o 


1 








-> 




r*H 










_J 




_J 


_J 


_J 








o 




1— 1 


to 















X 


X 


X 


X 


LU 


UJ 


UJ 


UJ 


UJ 


UJ 






X 


— 




D 








h— 




1 — 




1— 


1 — 


1 — 


> 


> 


> 


> 


> 


> 


X 






1— 1 




O 




r-H 




to 




lO 


(X 




















u 




to 


rv 




•— i • 




(O 












II 


II 


il 


X 


X 


X 


X 


X 


X 


< 




1 — 


r-H 


o 


> 




r> 




_J 






CO 








1 — 


1- 


b- 


H- 


1— 


1 — 


X 




lO 






UJ 




co 










o 


— 


— 


— 














1— 




UJ 


\— 


• 


X 






1— 


X 




X 


X 


(XJ 


r— I 


r\j 


II 


II 


II 


II 


II 


II 




h- 


h- 


to 


o 


X 


o 


II 


to 


h- 




1 — 


X 




«•— 


— 
















» — i 




►— ♦ 


UJ 






' 


►— t 








— 


1 — 1 


X 


X 




— 


— ‘ 


— * 


— ■* 


— * 


co 


O 


o 




• 


XT 


II 


“0 


_J 


II 


X' 


II 


II 


— 






co 


<r 


m 


vO 


r- 


00 


>— 


CO 


co 






u 










3 






1 — 






















co 


X 




< 


r-H 


(XJ 


X 


r“H 




(XJ 


1 — 1 


t0 






X 


X 


X 


X 


X 




UJ 


O 




1 — 




X 




co 


1— 




►—4 






*— 4 


u 


u 


u 


u 


u 


u 


u 


u 


I— 


1 — 


X 




w 


1 — 


CO 


co 


— 


no 




CO 


CO 


_J 


< 


< 


< 


< 


< 


< 


< 


< 


♦— * 




UJ 


II 






o 






o 




o 


o 




X 


X 


X 


X 


X 


X 


X 


X 


X 


o 


1— 




X 


. X 


X 


o 


X 


X 


o 


X 


X 


X 


\— 


1— 


1— 


1— 


1 — 


1— 


1— 


1- 


3 


o 


tO 


v: 


i — « 


*-• 


X 


o 


— 


X 


u 


X 


X 


1— 


1— 1 






















o 






1— 1 




<r 


in 


C\J 








(XJ 






















CO 






co 




co 


co 


co 








CO 






















co 






co 




co 


co 


co 









u u u 



u u u u u u 



102 



WRITE (4) TR LIST 
GO TO 3017 



SUBROUTINE 1000 PRODUCES A NEW TRACK HISTORY FOR EACH 
POSSIBLE CORRELATION IN THE TR LIST. 
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307 CONTINUE 
GO TO NEXT 



ALL TENTATIVE TRACKS ARE NOW ON UNIT 4. ELIMINATE DUPLICATES 
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DO 405 J = 1,1 
READ (4) TR LIST 
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PREPARE NEXT PERMANENT RECORD AND PRINT RESULTS. A KEY STOP IS 
PROVIDED HERE IN CASE A NEW TAPE IS DESIRED FOR PERMANENT RECORD 
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CALL EACH TRACK INDIVIDUALLY AND PRINT 
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900 END FILE A 
REWIND A 
REWIND 3 

9000 I F ( SENSE SWITCH 3) 90A, 9001 
90 A PAUSE 7000 












CO 

o 

o 



< 1 - 

o 

o 



o 

— r - 

o n- 
m oo m 

1 — I r*H O 

CO <\J 
' — I CD • 

• <\l I s - 
CO • lO 

vO 

• * 

h- I 



•—<•—<0 0 1 — 
O _J _J X 

O ^ LU 

O O II II 2 
-J 

o * — i *— < o 



2 e o 

O U. O O O 2 

O w j j O LU 

. c\j co <j- 

o o o 

o o o 



113 



FUNCTION DEGREES ( X » Y ) 
IF ( Y .EQ. 1.57079 ) 
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PROGRAM READOUT 
DIMENSION DATA! 100,3 ) 
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182 DAT A ( I » 2 ) = N + FRACTION 

183 X = ABSF ( DATA ( 1,3)) 
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APPENDIX IV 



FLOW CHARTS FOR THE PROBABILISTIC MODEL 

The flow charts in the following pages are intended to give the 
reader a general picture of the manner in which the logic of the 
probabilistic model was implemented in TRACK IV. Detailed flow charts 
may be obtained by contacting the Operations Analysis Department of 
the U. S. Naval Postgraduate School. 

The quantity 1^^ which appears in the flow charts is defined 
as the probability that a tentative track A.n established on sweep i 
is the true extension of track A in view of the possible extensions 
determined on sweep i+1. is computed as the maximum p b q over 

all possible extensions of track A.n multiplied by the track reliabil- 
ity of track A.n computed on sweep i. 

L is used to eliminate tentative extensions established on 
max 

sweep i which are duplicates of tentative extensions for other tracks 
established on sweep i. The duplicate with the highest 1^^ is re- 
tained. Ifrax is also used to determine which branch of a split track 
should retain the identity of the track and which should be redesig- 
nated. The branch with the highest 1^^. retains the former designation 
and all others are redesignated. 
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FLOWCHART SYMBOLS 




start 




stop 



computation or replacement 
procedure 




A 



DO loop 



termination of DO loop or junction 
of two or more branches of the 
program 




link between pages 




conditional GO TO 
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Note: All duplicate tentative 

tracks have been eliminated. 
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remaining TR LIST 



V 

For each 
correlation in 
TR LIST 



V 
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